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Abstract. Consider an N X N hermitian random matrix with independent 
entries, not necessarily Gaussian, a so called Wigner matrix. It has been 
conjectured that the local spacing distribution, i.e. the distribution of the 
distance between nearest neighbour eigenvalues in some part of the spectrum 
is, in the limit as N — > oo, the same as that of hermitian random matrices 
from GUE. We prove this conjecure for a certain subclass of hermitian Wigner 
matrices. 



1. Introduction and main results 

Consider a probability measure Pat on the space of all N x N hermitian matrices. 
We will be interested in the statistical properties of the spectrum as N becomes 
large, in particular in features that are insensitive to the details of the particular se- 
quence of probability measures we are considering. It is believed, on the basis of nu- 
merical simulations, that for many types of hermitian random matrix ensembles, i.e. 
choices of Pat, the local statistical properties of the eigenvalues are the same as for 
the Gaussian Unitary Ensemble (GUE), where d¥ N (M) = Z N X exp(-f Tr M 2 )dM. 

Here dM is Lebesgue measure on the space Hn ~ of all N x N hermitian ma- 
trices. The asymptotic eigenvalue density as N — > oo (density of states) is given 
by the Wigner semicircle law p{t) = j L -\/(4 — t 2 ) + . Let Pn{xx, ■ ■ ■ , xn) be the in- 
duced probability density on the eigenvalues. The semicircle law is the limit of the 
one-dimensional marginal density as N — > oo. The m - point correlation function 

(1.1) R^\xi,...,x m ) = , N ' . / p N (x)dx m+1 . . .dx N , 
is given by, @ ch. 5, @, 

(1.2) R£\x 1 ,...,x m ) = det(K N (x i ,x j ))% =1 , 
where the kernel Kj^ix, y) is given by 

(1 3) K N (x y) = PN{x)p N -i(y) ~ p N -i(x)p N (y) c _ N(x 2 +y 2 )M 

k n x - y 

Here pn(x) — k^x n + . . . are the normalized orthogonal polynomials with respect 
to the weight function exp(— Nx 2 /2) on R (rescaled Hermite polynomials). From 
these formulas, and Planchcrel-Rotach asymptotics for the Hermite polynomials it 
follows that 
(1.4) 

K» TW^^Hu + • • . ,« + ^h) = det (^ {U ~ ^ 

«■•>'>■! « ' , , | ...... i . ; / 



n^oo (Np(u)) m m v Np(uy Np{u)' V *i) 
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if p(u) > 0. It has been proved, Q, @, B, that this is also true in other invari- 
ant ensembles of the form dP N (M) = Z^exp(-NTrV(M))dM . The orthogonal 
polynomials in (1.3) are then replaced by polynomials orthogonal with respect to 
exp(— NV(x)) on K. That the ensemble is invariant means that the probability 
measure is invariant under the conjugation M — > U~ 1 MU 1 with a unitary matrix 
U. Sufficient control of the limit ( |l.4| ) for all m > 1, makes it possible to deter- 
mine the asymptotic spacing distribution, i.e. distances between nearest neighbour 
eigenvalues, see 0. More precisely, let {tjA be a sequence such that t^ — > oo but 
t N /N -> as N -> oo and define, 0, @, S N (s,x), s > 0, x G R N , to be the 
symmetric function, which for x± < ■ ■ ■ < x^ is defined by 

(1-5) S N (S, X ) = J-#{1 < i < AT - 1; - Xj < ^_ | Xi - „| < 

Given an hermitian matrix M let a?i(Af) < . . . xn(AI) be its eigenvalues; we write 
x(M) = (xi(M), . . .,x N (M)). Then it is proved in @ that 



(1.6) lim E N [S N (s,x(M))]= / p(cr)da, 

JV^oo J 

for a large class of invariant ensembles. Here p(cr) is the density of the (3 = 2 local 
spacing distribution, the Gaudin distribution, given by the probability density 

(1-7) p(s) = -^det(I-K) L i {0tS) , 

where K is the operator on L 2 (0, s) with kernel K(t, s) — sin7r(t — s)/Tr(t — s), the 
sine kernel, see [ p0[ . 

The aim of the present paper is to extend ( |1.4|) a nd ( |1.6|) to other, non-invariant 
ensembles. It is conjectured, see |^(| p. 9, that (|1.4|) and (|1.6| ) should hold also for 
so called Wigner matrices where the elements are independent but not necessarily 
Gaussian variables. In this case the probability measure is not invariant under 
conjugation by unitary matrices. For other results on Wigner matrices see for 
example |p|, fp[, p[, [p7| and pgj. In particular, in (26) the universality 

of the fluctuations of the largest eigenvalue is established. To be more precise, 
consider the complex random variables Wjk, 1 < j < k with independent laws 
Pjk = Pjl ® p jk, where p jj = s o- Let W p , a class of Wi gner ensembles, denote the 
class of all {Pjk}i<j<k which satisfy 

(1.8) J zdP. jk (z) = , J \z\ 2 dP. jk (z) = a 2 

for all 1 < j < k, and furthermore 



(1.9) sup / \z\ p dpj k (z) < oo. 

If Wkj — Wjk, W = (wjk)^ k=1 is an N x N hermitian Wigner matrix . 

Fix a > and let 4> a (t) = (7ra 2 ) -1 / 2 exp(— t 2 /a 2 ) be a Gaussian density function. 
Define Q* 1 = (f> a * P*/ , 1 < j < k, Q% = ^ * Pf v J > 1 and = S a . Then Q 
is also a Wigner ensemble and we let denote the subclass of W p obtained in this 
way. Note that although W% does not contain all Wigner ensembles it does contain 
cases where the distribution of the matrix elements have very different shapes, so in 
this sense it is rather broad, and proving universality in clearly shows that the 
universality is not restricted to the invariant ensembles. Another way to describe 
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this ensemble of random matrices is as follows. Let V be a GUE-matrix with the 
probability measure exp(-±TrV 2 )dV and let W be an N x N Wigner matrix 
with distribution P E W p , i.e. the law of Wjk is Pjk- We will assume that the 
variance a 2 — 1/4, which can always be achieved by rescaling. Then W + aV has 
the distribution Q, and we write 

M = -$=(W + aV). 



We can think of this in terms of Dyson's Brownian motion model, M , W + aV is 
obtained from W by letting the matrix elements execute a Brownian motion for a 
time a 2 , see sect. 2. If P £ W p and W is is an N x N hermitian matrix we let 
PW denote the distribution of H = W/\/N = (hjk), i.e. 

dPW(H) = 11 dP jk (^Nh jk ). 

!<3<k<N 

The matrix M has the distribution Q^ N \ which is given by 

(1.10) dQW(M) = 2~ N ' 2 (JLy /2 ^J e-^( M - H ?dPW(H)jdM, 

and this is the measure we will study. The asymptotic distribution of the eigenvalues 
xi, . . . , xm of M is the semicircle law 

2 



(1.11) p(u)= 7r(i + 4a2) V(l + 4« 2 - M 2 ) + . 

The following proposition will be proved in sect. 2 using an argument from 0], [H- 

Proposition 1.1. The symmetrized eigenvalue measure on induced by Q^ N ^ 
has a density 



(1.12) Pn(x)= p N (x;y(H))dp( N \H) 

where 

( , 13) »<*v)- (£)"'' 

and An(x) = Yli<i<j< n ( x * ~ x i) * s ^ e Vandermonde determinant. 

The main result of the present paper is that for Wigner ensembles from W p we 



can prove (1.4) and (1.6), and thus extend the universality to a rather broad class 



of Wigner matrices. 



Theorem 1.2. Fix a > and assume that \ u\ < y/l/2 + 2a 2 . Let R JJ {x\ 



be the correlation junctions, defined by (1.1), of the eigenvalue measure pn, 



for QW, ( |jJ6p . Let f € Lf(M. m ), the set of all L°° functions on R m with compact 
support, and set for x £ M. N 

(Sf)(x)= /fe!-"^*™). 
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where the sum is over all distinct indices from {1, . . . , N}. If Q S W% with p > 
2(m + 2), then 

(1.14) lim / {Sf){Np(u)(x 1 {M)-u),...Np(u)(x N {M)-u))dQ {N ' ) {M) 

TV— >oo 



n^oc J Rm (Np(u)) m Np(u) Np(u) 



f( tl , ...,t m ) det C^ y )% =1 <rt. 
n(U - tj) 



The condition on u is made just to simplify the saddle-point argument in sect. 
3; the result should hold for any u with p(u) > 0. 

We can also prove that the spacing distribution is the same as for GUE. 

Theorem 1.3. Fix any a > and assume that Q £ Wf +e , e > 0. Let Sn(s,x) be 
defined by Then, for any s > 0, 

(1.15) lim f S N (s,x{M))dQ( N) (M)= [ p(a)da, 



Hi 



where p{s) is given by ( |_?.'/| j. 

The theorems will be proved in sect. 4 after the preparatory work in sect. 2 and 



3. 



2. The correlation functions 



We will start by proving Proposition 1.1 using the Harish-Chandra/Itzykson- 
Zuber formula following Q , Q . After that we will give a formula for the correlation 
functions of pN{x;y), which is very close to the formula in H, but our derivation 
will be different. A central role will be played by non-intersecting one-dimensional 
Brownian motions and we will use the formulas of Karlin and McGregor. Also we 
will discuss the relation to Dyson's Brownian motion model. This connection can 
be found in [^) and we will only give an outline. 

Proof. Let F(x) be a continuous symmetric function on M . By Fubini's theorem 
(2.1) 

F{x{M))dQ^{M)=c^ f (f F(x(M))e-^ Tl{M - H)2 dM\dp( N \H) 

with e)y = 2~ JV / 2 (A ~/ira 2 ) N ' '/ 2 . In the right hand side of (2.1) we make the substi- 
tution M = U^RU, with U € U(N) and R S H N , and then integrate over U(N). 
If we use Fubini's theorem again, we obtain 



F(x(R)) / e-^ Tl ^ lRU - H ^dU dR dP( N \H). 



Hi 





Here we have also used the fact that dM = dR. The integral over U (N) can now 
be evaluated using the Harish-Chandra/Itzykson-Zuber formula, , |l3[| , see also 
|20| A. 5. We obtain the integral 

c « c ( 2) f f f BY^ - Art.( f .- N <.*i-v*) a l*> 



I (Ik M dfit(e-^-»*) 3 / a » a )y t=ldfl ) dP^m 

Jh n A N {x)A N (y) J - J 



where y±, . . . , i/n are the eigenvalues of H and = [a 2 /N) N ^ N ~^/ 2 rijLi 3-- The 
integrand in the middle integral depends only on the eigenvalues x of R and hence 
we can integrate out the other degrees of freedom in the standard way, ]2(| ch. 3, 
and obtain, after using Fubini's theorem, 



F{x(M))dQ {N) (M) 



(2.2) 



J1)J2) „(3) 
N N L jV 



Hn \JR N 



F(x)^l det(e- N ^-y^/ 2a2 )"d N x) dP^(H) 



with cf = Tt N{N - 1)/2 Y[" =1 {j\)- 1 . We see that c^c^cf = {N/2^a 2 ) N / 2 and 
since (|2.2|) holds for arbitrary bounded, continuous, symmetric F{x) we have proved 



that the symmetrized eigenvalue measure is given by 1.12. This proves proposition 

m □ 

Let pt(x,y) be the transition probability of a Markov process X{t) on R with 
continuous paths. Consider N independent copies of the process (Xi(t), . . . , Xjv(t)) 
and assume that this is a strong Markov process in M. N . Suppose that the particles 
start at positions yi < • • • < yN at time 0. The probability density that they are 
at positions x% < ■ ■ ■ < xn at time S given that their paths have not intersected 
anytime during the time interval [0, S] is, by a theorem of Karlin and McGregor, 



|16|, 



det(P s ( %) * fc ))£ fe=1 . 

Hence, the conditional probability density that the particles are at positions y\ < 
■ ■ ■ < yN at time 0, at positions x\ < ■ ■ ■ < xn at time S, at positions z\ < ■ ■ ■ < z/v 
at time S + T, given that their paths have not intersected in the time interval 
[0, S + T] is 

(2.3) q s . T (x;y;z) = ■^-det{P s (y J ,x k ))f k=1 det{P T (x j ,z k ))f k=1 , 



N 



where 



Z N = det(P s (y J ,a; fe )) Jife=1 det(PT(a:j,z fe )) ifc=1 d x; 

J x\<---<x N 



we assume that Zn > 0. Note that the expression (2^3) is a symmetric function of 
xi, . . . , xn, so we can regard it as a probability measure on K. . Our next lemma 



shows that we can obtain p^(x;y) defined by (1.12) as a limit of the measure in 



(2_a) 



Lemma 2.1. Let Zj = j — 1, 1 < j < N and let pt(x,y) — (2nt)~ 1 / 2 exp((a; — 
y) 2 /2t) be the transition probability for Brownian motion. Then, for any x € M. N 
and yi < . . . ,y N , 

(2.4) Umwtowz) = ^WfiZU ^'^^ 1 = qSiX;V) - 
Note that p N (x;y) = q a 2 /N (x;y). 
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Proof. Write 

dBt{Ps( yj ,x k ))? tk=1 det(P T (x i! 

1 N «=!+*? 

(2-5) = (^4^ det(e " (Xj "" )2/2g)3Vl n e "^ det^-/-)^!- 

Note that Zjv is the conditional probability density of going from y\ < . . . ,yjv to 
z\ < . . . , Zn without collosions, i.e. 

Z N = det(p s+T (y v z k ))f k=1 

1 N y 2 +z 2 

(o 6) - - TT e^^srt) det(e y i Zk/2{s+T) ) N , 



(2n) N / 2 (S + T) N / 2 l \ 



j,k=i- 



Now, since Zj = j — 1, we have two Vandermonde determinants in (|2.5|) and (2.6) 



If we evaluate these, take the quotient between (|2.5|) and (|2.6D and then take the 



limit T — > oo, we obtain the right hand side of (2.4). □ 



Proposition 1.1 and lemma 2.1 establish a link between the eigenvalue distribu- 
tion of M = (W + aV)/^/N and the non-intersecting Brownian paths. If we set 



5 = a /N, then the right hand side of (2.4) and (1.13) are identical; yi < ■ ■ ■ < un 



are the eigenvalues of H = W/v N. This relation can also be seen in another way, 
which we will now outline. Let X(t) = (xj k (t))^ k=1 be an N x TV Hermitian matrix, 
where Kexj k (t), Im Xj k (t), j < k are independent Brownian motions with variance 
(1 + 5jk)/2. Assume that X(0) — H is distributed according to P( N \ Then the 
distribution of X(o? JN) is the same as that of M = (W + aV)/^f~N. Following 
Dyson, |^|, see also |2J|, it is possible to derive a stochastic differential equation for 
the eigenvalues Ai(i), . . . , Ajv(£) of X(t), 

(2.7) dX r = dBi + V - l — dt, 

where Bi are independent standard Brownian motions on R, and with the intial 
conditions Ai(0) = j/j, 1 < i < N. We can also consider the problem of non- 
intersecting Brownian motions in a different way than that of Karlin and McGregor. 
Namely, let K = {x E R w ; x\ < ■ ■ ■ < x^} and consider Brownian motion in M. N 
starting at y S K and conditioned to remain in K forever. As proved in ]^], see 
also [O] , p4| , if Xi are the components of the A-dimensional conditioned Brownian 



motion they satisfy the stochastic differential equation (2.7) with the same initial 



conditions. This gives another way to obtain (1.13) without using the Harish- 
Chandra, Itzykson/Zuber formula. Actually, we can turn the argument around 
and give a proof of this formula. 

We turn now to the computation of the correlation functions of the right hand 



side of (|2.4|), but we start more generally with (|2.3| ). This can be analyzed using the 
techniques of |2£|, compare the analysis of the Schur measure, ^l), in ji~4|], and see 
also [[[I]. For completeness, let us outline the result we need from ^9). Let (f2,/i) 
be a measure space. Assume that <f>j,ipj G L 2 (Q, /i), 1 < j < N, and / S L°°(Q, y). 
Set 

N 

Z N [f] = W[ / det^-Orfc))^ det(^(x k ))» k=1 /(a^d/ifo) 
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and 



A = 



(j) j (x)tp j (x)dfi(x] 



N 



3,k=l 



Proposition 2.2. (|29[]J Assume that Zn[1] =/= 0. Then A is invertible and we can 
define 

N 

(2.8) K N (t,s) = MWA-^kjM*)- 

3,k=l 

Then, for any g € L°°(Q, ji), 

Z N [l+g] 



(2 ' 9) z sm 

If we define a density on n, N by 

1 



(2.10) 



u N {x) 



N\Z N [l] ~" v ' 
then it has the correlation functions 



= det(I + K N g) L 2 {n) . 



det(^(x fc ))^ fc=1 det(V'j(a; fc ))f fc=1 , 



(2.11) 



{N - M)\ J UN 



upf(x)dx 



m+l 



,dx N = det(K N (x l ,x j )) ij=1 



Proof. We will indicate the main steps in the proof of (2.E) of which (2.11) is a 
consequence, see 29 . Set 

N 



D 



(f)j{x)tpj{x)g(x)d^{x) 



j,k=l 



Then, by the formula 



Z N [f}= det [ / cj )j (x)^ J (x)f(x)d^x) 



N 



3,k=l 



which goes all the way back to (TJ , and which is not difficult to prove by expanding 
the determinants, we see that det A — Zn[1] 0, so A is invertible and 

Z N [l+g] _ dct{A + B) 



(2.12) 
Now, 



Z N [1] 



det A 



N 



det(J + A~ B). 



(A- x B) jk = jj k (x) l^A^yMxM^j d»(x), 

and we define T : C N — > L 2 (Q,dfj,) and S : L 2 (Q,dfj,) -> C N by the kernels 
T{x, k) = tpk(x) and S(j,x) — J2e'=i(A~ 1 )je4 l e(Xg(x) ■ Then, by (2.12) and a deter- 
minant identity, 

z N [i + g. 



z N [i] 



det(I + ST) C n = det(J + TS) L 2 
det (J + K N g) L 2 {n ^ )} 



with Kn given by (2.8). Note that -Kjv<7, which means first multiplication by g 
and then application of the operator on L 2 (fl, d[i) with kernel Kn, is a finite rank 
operator. □ 



8 



K. JOHANSSON 



Observe now that if we take O = R, d[i(x) = dx, cj)j(x) — Pt{x, Zj) and ipj(x) = 
ps(yj,x), then (2.3) is a probability density of the form (2.10) and we can apply 
the proposition. Note that 

(A) jk -- 



p T {x, Zj)p s (yk, x)dx = Ps+riVk.Zj). 



The kernel which gives the correlation functions is 



JV 



N 



K^ T {u,v) = ^2p s (yk,v) ( ^2(A ^jkPT^Zj) 



fc=i 



Let Ak(v) be the matrix we obtain from A by replacing column k by (px{v, z\) . . . pr{v, Zn)) T ■ 
Then, by Kramers' rule, 



(2.13) 



Kf/ T (u,v) = ^2ps(Vk,v) 



k=l 



detA k (v) 
detA 



This formula and proposition 2.2 is the basis for the next proposition. The result 
is closely related to the result derived in M by different methods. 



Proposition 2.3. The correlation functions for qs(x;y) defined by (2^) are given 
by 



(2.14) 

where 
(2.15) 



;y) = 



(n -my. y R «- 



q s {x;y)dx. 



m+l 



, dx 



N 



KN( u ,v;y) 



{v - u)S(2iri) 2 
w + z — v — S 



dz I dw(l - e {v - u)z/s ) 

y.i 



(w -2vw- z +2uz) /2S 



(w - yj)(z - Vj) 



Here 7 is the union of the curves t —* —t + iu>, t g R and t —*t — iuj, t S R with a 
fixed uj > 0, and T : R 3 t — ► it. 

Proof. We have to show that with p t (u,v) — (27 rt)~ 1 / 2 exp(— (u — v) 2 /2t) and 
Zj = j — 1 the limit of the right hand side of (|2.13| ) as T — ► 00 can be written as 
( [2.15| ). The result then follows from lemma 2A, proposition and the dominated 
convergence theorem. We see that 



(2.16) 



dct .4 = 



(2tt(S + T)) N /'- 



N 

n 



* 2 +y 2 
3 3 



g 2(S + T) 



l<i<j<N 



by the formula for a Vandermonde determinant. Let be the curve t — > t + iM , 
t e R, M fixed. Then 



p T (zj,v) 



1 



V2^T 



_g 2(S + T) 2T 



1 



T 2 +2j(y+ir^/ 2 T(S + T))^ T 







Hence, 

' JV 



x e 2T / e T 2 det Afe(u)dT, 



where Ak(v) is the matrix we get from (exp( ifflj, )) jfe— i by replacing column A; 

by (exp(zj(^ + ir^/ 2T(J+T))))j=i- Since Zj = j - 1 we have a Vandermonde 
determinant and we obtain 



1 /' x 2 T r 



(2.17) xe^-p / e~~ II (e^+r - e^+ r )dr, 

V 21 l<i<i<iV 

where should be replaced by (5 + T)(y + * T ^/ 2T(J+t))- Take the quotient of 
( fjj ) and (|2.17D and let T -> oo. This gives 

hm det^r,)! /• e _ 4 „ ft, + fr/Sr - gJ V 
t^oo dot A V2^Jr* M j£\ Vk-Vi J 

Choose M so that v—ySM — L, where L is given, and make the change of variables 
w = v + i^/Sr. Then 

det Ak(v) 1 f (™-v) 2 tt / w — Vi \ , 

lim — - h ) ' = - / e~^~ TT — dw, 

t^oo det A iV2^SJr L L ,\ \Vk - Vj J 

where r L : t -> L + it, i e E. Thus, using (|2.13|) , 

Let 7 be a curve surrounding y\ , . . . , j/jv and choose L so large that 7 and T do not 
intersect. The residue theorem gives 

N2/oo AT JV 



— / e ~ (z " M) 725 TT = y e-(— ) 2 /^ TT f 

1 j=i yJ fe=l j^k vy yj 

for all aieTi. Thus, 

(2.18) K s N {u,v-y) = ^-(dz( dwe ^ 2 -^ w -^+2 UZ ) _J_ TT 

(2 7 rj) 2 S' y 7 7 ri w - 2 / = i z ~ % 



In (2.18) we make the change of variables z bz, w bw with tgK close to 1. 
This will modify the contours but we can use Cauchy's theorem to deform back to 
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7 and T^. Now, take the derivative with respect to b and then put 6=1. This 
gives the equation 



= K§(u,v;y) + 



v 2 -u 2 

e 2S 
(2ni) 2 S 2 



1 



7 ^ «» - * 

JV 



—2vw — z +2uz) 



w 2 ~ Z 2 + UZ — VW + S 



-\w- yj z-yj 



N 

n 

3=1 



This can be written 
d 



v 2_ u 2 

e as 



7 "/rt 



AT 



w + z 



*Jj 



{ (u> - Vj)(z - Vj) 



N 



(w 2 -2vw-z 2 )/2S e uz/S TT W ~ V] 

fA*-vi 



and integration of this formula gives (2.15). In this last formula we can choose L 
arbitrarily and take 7 to be the curve in the proposition by using Cauchy's formula, 
This completes the proof. □ 

We now take S = a 2 /N and set 
(2.19) K N (u, v; y) = e^^+^^V v; y), 

where w is a const ant that will be specified later. Note that we can replace ^ N 
with /Cat in ( 2.14 ) without changing the correlation functions, so we can just as 
well work with JCn- Set 



3=1 

if 2 N 



h(z, w) 



a u(u—v) 



Np(u)(v — u) 



(w-yj){z-yj) 



e j?( u ~ v ) w ( e N(u-v)w/a 2 _ e N(u-v)(w-z)/a 2 ^ 



dz f dw 



\h{z, w)g N (z, toJe^M")-/* W). 



so that 

(2.20) }C N (u,v;y)=Np(u) 

These are the formulas we will use in the asymptotic analysis. A straightforward 
computation shows that 



(2.21) 



g N [z,w) = -f N (z) + . 

z z — w 

3. ASYMPTOTICS 



The eigenvalues yx, . . . , yN of the Wigner matrix H converge to the semicircle 



law 
(3.1) 



a(t) = ~^/^^^ \t\ < 1. 
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In order to be able to perform the saddle point analysis of ( 2.20| ) we need uniform 
control of the convergence of Jn{z) to its limit 



(3.2) 



i r 1 

f(z) = ^(z 2 -2uz) + J log 



- t)a(t)dt. 



In order to show this we must start with some probability estimates. Write = 
{z e C; |Rez| < R,rj < \Jmz\ < R}. 

Lemma 3.1. Let F € L°°(R N ) be symmetric and let r\ > and R > be given. 
Assume that P € W p , p > 4 and < £ < min(i — |, i). Then, there is a 

probability measure pW on Hn such that 



(3.3) 

and 

(3.4) 



F{x(H))dP w (H) 



Ht 



F(x{H))dP( N) {H) 



Ht 



< jV 2 - p ^- fl ||F| 



sup 



i r 1 

-Trlog(z-H)- J log 



(z - t)a(t)dt 



< ON'* 



a.s. with respect to P( N ^ . 



Proof. Given P £ W p we introduce a cut-off L > and define a new probability 
measure Pl S W p by 

1 



where c?l is a normalization constant. Note that Pl jk is supported in K 



[-L,Lf. Set > =Ui< j <k< N dL, j k. Then, 
/ F{x(H))dP (N \H) - [ F{x{H))dP { L N) {H) 



<\\F\Ul-d^)(l 



(3.5) 

for some constant C. The last estimate follows from 



CN 2 



(3.6) 



1-d 



(N) _ 



P[some \W jk \ >L]<N 2 sup 



E[\W jk \v] < CiV' 



1<3<* 



LP 



by ( |l.9| ). Set = f^,, n -^Z 2 and note that #Dn < CN 2 for some constant C 
that only depend on R, n. For a given function / set 

A N (f;S) = {HeH N ; \±Tr (f(H)) - J J(t)d*(t)\ < 5}, 



where a(t) is the semicircle law (p.5\). Set 



(3.7) 



where f z (t) = log(z — t) (principal branch). To estimate the probability of An (8) 
under pW we will use a result of Guionnet and Zcitouni, fliol. Let 



\f\c = sup 

t,s£X 



\x-y\ 
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and ||/||£ = H/Hoo + \f\c- Then, by [|10|, corollary 1.6a), and the discussion before 
this corollary, given e > 0, there are positive constants Co(e), C\ and C2 such that 
if we write 



(3.8) 
then 

(3.9) P 



S^N) = C^lflcN- 1 + C^WfWcN- 1 ^ 



(AT) 
L 



\±Tr f(H) - £j(t)a(t)dt\ 



> 8 



< 4exp 



C 2 N 2 



(S-S 1 (N)Y 



for any S > Si(N). Since under P£ all |£/jfc| < ^/2{L/\fN), the spectral radius 
is < 2L. Thus, the left hand side of (3J3) is unchanged if we replace f = f z with 
/ = (t), where / z L (i) = log(z - t) if \t\ < 2L, f^(t) = log(z - 2L) if t > 2L and 
/f'(i) = log(z + 2L) if t < —2L. Now, (t) is Lipschitz and there is a constant 
C 3 , independent of L, such that \f^{t)\ c < C 3 and WfH^Wc < C 3 (l + logL) for 
all z G n Z) ,. Take L = L N = TV 1 / 2 "? and e = 1/6 in Q. Then Si(N) < CN~ 2 ^ 
and if we choose S = iV« in (O) we obtain 



(3.10) 



P, 



(iV) 



|iTr/,(iJ) - I /,(t)o-(t)dfcj > AT« 



< ciexp(-c 2 iV 2f ) 



for some positive constants c\ , c %. I f we use ( 3.10] ) we see that the probability of 
the complement of the event in ( |3.7|) can be estimated as 



(3.11) 
Set 



P { L N) [A N {N-ty] < CN 2 e- c * N2e 



dPW{H) = (P i L N N ) [A N (N-i)])-' XAN(N - () (H)dP J 



(JV) 



Note that N 2 /L P N = jv2-p(i/2-« j so combining (p^), ( p7| ) and ( pd| ) we obtain 
the estimate (3.3). From the definition of An (6) we see that (3.4) holds for z G Dn, 
but then a straightforward approximation argument extends it to all z € £Ir, v . This 
completes the proof of lemma 3.1. □ 

We now come to the central asymptotic result. 

Lemma 3.2. Let fi_R )?7 be as above, let £ G (0, 1/2] and let K be a compact subset 
ofM.. Also let un be a sequence such that un — * u as N 00. Furthermore, let 
Yr,t) be the set of all y £ M. N such that 



(3.12) 



sup 



1 N 



z-Vo)- I lo g( z - t)a(t)dt 



< CN~ 6 



for some constant C and all N > 1, where a(t) is given by (3.1). Then, we can 
find Rq > 0, 770 > and a constant C such that for all y £ YRo,r) > T € K> 
\u\ < y/l/2 + 2a 2 and N > 1, 

1 r 
( 3 - 13 ) at 1 -, ^n(un,un + — , s \y)- 



Slll 7TT 



7Vp(u) v ' Np(u) 
where p(u) is given by ( | 1 



< C(\u-u N \ + N-t), 
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Proof. It follows from the formula ( 2.20| ) that 



(3.14) 
1 



Np{u) v Np(u) ' L 2ni 

where gpf(z,w) is given by ( |2.21 ) 



Jn(z) = (z 2 ~ 2u N z) 



—h(z,w)g N (z,w)e N ^^'^\ 
2ni 



1 N 

3=1 



and 



(^ e -Tw/a 2 p(u 



-t(w — z) I a p(u) 



h(z, w) = - 

We have taken u> — ujqJN p{u), where ui is given by (3.23) below. The integral 
in (3.14) will be analyzed using a saddle point argument. It follows from (3.12) 
and Cauchy's integral formula that there is a constant C such that for all TV > 1, 
TEK,ye Y R/2av and |u| < ^1/2 + 2a 2 , 

(3.15) \f N (z) - f(z)\ < C(N-t + \u- u N \) 

\fMz)-f"(z)\<CN-e. 

A computation shows that, 



f(z) = -Az -u) + 2(z - Vz 2 - 1). 



(w + l/w)/2 with inverse S 1 (z) 



z + \/ z 2 — 1, where \/ z 2 — 1 



S et S(w ) 

\/z — \\Jz + 1 (principal argument). The function S maps {\w\ > 1} to C \ [— 1, 1] 
and \w\ = 1 is mapped to [—1,1]. Note that 



f'(S(w)) = 



2a 2 



(2 



1 1 v_ 

2a 2 w a 2 



Write u = yl + 4a 2 cos 9 C , where 9 C S [0, ir]. Our assumption on u means that 
I cos6» c | < 1/2. Note that f(S(w)) = has the solutions wf = y/l + Aa 2 exp(±i0 c ). 
Hence the critical points for / are zf = S(wf). 

We will now define some contours that we will use. Pick S > (small), see 
below. Set, for some e > (small), jf(t) = S{\/1 + 4a 2 e iS - t), -oo < t < 0, 
7 +(t) = 5(Vl + 4aV t ), 5 < t < 6 c -e, y£{t) = S{Vl + 4aV t ) , 9 c -e< t<6 c + e, 
7+(i) = S(Vl + 4a 2 e it ), 9 C + e < t < tt - S and 7 +(f) = 5(Vl + Aa 2 e 1 ^-^ - t), 
< t < oo. Also, set 7j"(i) = Jj'it), 1 < j < 5. Then, we can take 7 = 



£,-=i(7j - 7, ) 
ImS(A 



7" 



+1 - 



in ( j3.14| ). Let i S (l/\/l + 4a 2 , 1) be such that 
= 77, and write a = ReiS , (ioW+). Set, for some e > (small), r^~(t) = 
a + it, < t < rj, T+(t) = S(tw+), t < t < 1 - e, T%(t) = S(tw+), 1-e < t < 1 + e 
and T+(t) = S(tw+), 1 + e < t. Also, set iy(i) = r+(i), 1 < j < 4. We can then 

taker = £ 3 l 1 (r+-ry) = r 



(3.16) 
where b, d G { 



T+ - r- in ( |3.14| ). Set 

L>tf(T;y) = N [ ^[ pLh(z,w)g N (z,w)e 
J b 2m J T d 2m 



-} and write L 



N 



N 



L 



N 



J N 



N(f N (w)-f N (z)) 



J N 
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Claim 3.3. We can choose Rq > 0, r)o > ande,5 > 0, so thatj 3 +-f 3 +r^ + T 3 
lies in a neighbourhood of z^r which is included in fi_R /2.2tj an d f or a ^ N > 1, 
t £ K,y G Y r/2 ,2n and M - V 1 / 2 + 2a ' 2 > 
1 



(3.17) 

with c > 



Np(u) 



ICn(un,u n 



Np(u) 



y) - L N (r;y) 



< Ce 



-cN 



The cl aim will be proved b elow. We will now use the claim to finish the proof 

N 



of lemma |3.2|. It follows from ( |3.15| ) that there are critical points z^ = S(w^) for 
/n(z) such that 

(3.18) \z±-z±\<C(N-e + \u-u N \). 

We can deform 7^ (T 3 ) into contours 7^ (r^) such that the endpoints are un- 
changed, 7^(0) = r^(0) = z± and 7^ (r± ) have C^-distance < C(iV~ s + \u-u N \) 
to 7* (Tf ). We can also asume that these contours are chosen so that 7#(t) = 
S{w±e ±lt ) and Y%{t) = S{w%{\ + 1)) for \t\ < e. 

We can now proceed in the standard way with a local saddle point argument in 
( [3.16| ) and prove that there is a constant C such that 



(3.19) 

L b *( T ;y) - h(z b N ,z d N )g N (z b N ,z d N ) 



2tt 



(7^)'(0)(r^) , (0)e JV W«( z «)-/«( 2: «)) 



(2 ™ )2 v / ^(^)(^) , (o) 2 v -/^(4)n)'(o) 2 



< 



c_ 

N 



-N 



for all N > 1, r G K , y G Y Ra ^ and \u\ < y/l/2 + 2a 2 . Note that 
and f N (z^f) - Jn{z n ) is purely imaginary. Now, (7^)'(0) = biS'(w b N ), (r^)'(0) 
w^S' (w 1 ^) and a computation shows that 

/^(4)(7^)'(o) 2 = -f'M)(r b N )'(o) 2 = -a{z b N )s'{w b N f{w%f, 



-N 



which has a positive real part by (3.15) and the fact that f"(z b )S'(w b ) 2 (w") 2 
has a positive real part. From ( [2.21 ) we see that Qn(^ni Z%) = if b ^= d and 



g N (z b N ,z b N ) = lt follows that 



M4r>4)(7&)'(0)(lV(°) 



V/^(4)(7^)'(o) 2 v /-/^(4)n)'(o) 2 



-bi. 



Also, from ( ft.l§l ) it follows that 
thus (|3.1S|) yields 



(3.20) 

if b =/= d and 
(3.21) 



b N ,z b N )-h(z b ,z b )\ <C(N-t + \u-u N \), and 



^(r;y)l< 



bh(z b ,z b c ) 
2m 



N 



< C(N-Z + |u — t*jv|)- 
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Combining fU6j ), ( |3.20[ ) and we obtain 



(3.22) 
Now, 



Ln(t;v) + 



2m 



< C(N~ S + \u-u N \). 



and a computation shows that 

z± 1 + 2a 2 



(3.23) 
Thus Q3.22D becomes 



a 2 p(u) 



2a 2 



sm tit 



COt 9 C ± 7TZ = Ll>0 ± 7TI . 



< C(JV-« + lu-ujvl). 



If we combine this estimate with (3.17) we see that the lemma is proved. 



□ 



It remains to prove claim 3.3. 
Proof. Let 7^ = if an d — J2j^3 ^f- We have to estimate 

jf = iV / \dz\j |d W ||/ l (z, W )|| 3w (z, W )|e^ Rc ^^^"^»- ArRo ^^-^^», 



7; 



and 



I b 2 d = N \dz\ / \dw\\h(z,w)\\g N (z,w)\e NKc ^^-^^~ NRc ^^-f-^\ 

where b,d G { + ,—}. Note that /jv(^) — /jv(^g — ) is purely imaginary. We will 
concentrate on since the other cases are similar. 
Using the inequality 



w-yj 



Vi 



<l + C(\w\ + \z\) 



it is not difficult to see that there are constants C\ and C2 such that 
(3.24) 

\h(z,w)\\g N (z,w)\e NRc{fNiw) ~ fN{z)) < C 1 £'C2JV(| Z | + |, i >|)+JV(Ro( lU 2 -2t l «,)-Ro (z 2 -2« Z ))/2a 2 



for all y G R N , r G K and \u\ < \A/ 2 + 2 ° 2 - Note that |Im z| > c> for all 2 G 7. 
(The constant c depends on the S in the definition of 7, but as we will see below 
6 depends only on the parameter a in the problem.) From the estimate (3.24) it 
follows that by picking R = Rq sufficiently large, the contribution to + from z 
and/or w outside ^Ir 0i q is < e~ N . Thus we can assume that z,w G Qr ,o- Next, 
we will derive the other estimates we will need to prove the claim. 
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Assume that z £ ^r ,jj and w £ Tj. Then, 
\g N (z,w)e N f"^\ 



< 



< 



D NRc (w 2 -2uw)/2a 2 



< (7g A '[ Re /»( a + i l)+ Rc (w' A -2uw)~((a+ir)) 2 -2u(u+ir)))]/2a 2 



If we use ( 3.12 ) and the definition of /n we obtain 

\g N {z,w)e N ^ N< -^- fN ^\ 

(3.25) < C e cAr ( JV_5 + l«-"ivl)+^ 2 /2a 2 +A r Ro(/(a+ir,)-/( Z +))/2a 2 

for z S ^K 0l r; and to G rj~. 

We will now compute how Re/(z) changes along 7. Assume that 9 C > 0, the 
other case is analogous. Consider 7(0) = S(yl + 4a 2 e l9 ), <5 < 6* < 7r — <5. A compu- 
tation, using the fact that f'(j(d c )) = gives Re^/(7(0)) = 
cos 8) . From this we see that there is a constant Co > such that 

\2 



(3.26) 



Re (/(Vl + 4aV e ) - /(*+)) > c (# - C ) 2 



Next, consider 7^) = S(Vl + Aa 2 e lS -t),t< 0. If we write u s = \/l + 4a 2 e iS , 
then 



dt 



4a z ujg — t 



1 



(w 4 - i) 2 



Set — < = s{t)e lS ( l \ A computation shows that 



(3.27) Re -|/(7i (*)) = -— — L= 
at 4a 2 Vl + 4a 2 



1 + 4a - 



1 



■cos20(t) 



(s(t) + — -)cos6(t) - 2 cos < 
s(i) 

-i-sin2^)( s (t)--^)sin^) 



Note that sin#(i) = s(t) 1 Vl + 4a 2 sin 8. It follows that the right hand side of 
fl3.27j ) equals 

(3.28) 



1 



-2 



4a 2 Vl + 4a 2 
(1 + 4a 2 ) sin 2 8 



(1 + 4a 2 ) sin 2 8 

2 7-71 ) 



«(0 4 



cos0(i) - 2(1 + 4a 2 



1 



„(1 + 4a 2 ) sin 2 (5 N 
2- r^; cos* 



S {ty 



and this is 

1 



< 



F (l + 4a 2 



1 



4aVl + 4a 2 ^ ' "' s{t) 2 
since s(t) > 1. Choose <5 < # c /4 so that 

1 + 4a 2 



(s(t) + 4t) cos " 2(1 + 1 » sin 2 (5) cos C 
s(i) 2a^ 



(1 



sin 2 8) cos # c < cos — . 
2a 2 ' ~ 2 
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Since s(f) + l/s(t) > 2 and 0(f) < 6 we see that there is a constant cq > such that 



(3.29) 



Re-/( 7 iW) <-co- 



For 75(f) = + 4:a 2 e i( - n -^ - f), f > 0, we still have the formula (|3.28|) with 

7i(f) replaced by 75 (f) and, since 7r — S < 0(f) < n, we see that the right hand side 
is 



(3.30) 



Vl + 4a 2 s(f) 



and consequently there is a constant cq > such that 



(3.31) 



Re-/( 75 (t))>co. 



Consider now how Re /(if) changes along T + . Set r(f) = S(tw^), f > fo- A 
computation gives 



' " t [1 + f(l + 4a 2 ) - (f 2 (l + 4a 2 ) + -j-)cos20 c 



Now, since \u\ < y/l/2 + 2a 2 , it follows that cos 20 c < and thus 



(3.32) 



Re j t f{S{tw c )) > ^(1 + f(l + 4a 2 )) if f < f < 1 
Re j t f(S(tw c )) <±Z±(l +t (l + 4a 2 )) if f > 1. 



The first of these estimates can be used to show that if we pick r\ = tjq sufficiently 
small, then 

7 ] 2 + R C (f(a + i V )-f(zt))<-c 
for some positive cq. If we use this in ( |3~25| ) we obtain 
(3.33) \g N (z, w)e N ( fN M ~ fN \ < Ce~4> N 

for some positive c' . We can now use ( 3.26 ), ( 3.29| ), ( |3.31 ), ( 3.32| ) and ( 3.33 ) to 
estimate + and see that it is < Ce~ cN for some positive c. □ 



4. Proof of the theorems 



We start with the proof of theorem 1.2 



Proof. By proposition 1.1 and Fubini's theorem the integral in the left hand side 
of (1.14) can be written 

PN (x, y(H))(Sf)(Np(u)( Xl - u), . . . , Np(u)(x N - u))d N x] dP^(H) 



>H N \JR N 

Note that H^/)!^ < A^H/Hoo- Since p N (x, •) is a probability density on R N we 
can use lemma 3.1 to replace the expression in (4.1) by 

(4.2) 



PN (x, y(H))(Sf)(Np(u)( Xl -«),..., Np(u)(x N - u))d N x dP^(H) 
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with an error < CN m \\f\\ 00 N 2 ~ p ^ 1 / 2 ~^ = o(l), since p > 2(m + 2), provid ed w e 
choose £ small enough. Now, since p^(x,-) is symmetric it follows from (1.13), 
(EO) , (|2.14|) and (|2 . 1 9j) that the expression in dO) can be written 



(4.3) 



H, 



f(tl, • • ■ , i m ) 



x dct( 



Np(u) 



JC(u 



Np{u) 



Np{u) 



y{H)W 1=1 d m tdP^ N \H). 



Since / has compact support and we know that (3.4) holds a.s. [P*^] it follows 
from lemma |3.2| , with un = u + ti/Np(u), r = tj — U, that 

sin ir (U — tj) 



JC(u 



Np(u 



Np(u) 



;y(H)) 



< CN~t, 



for a. a. [P^] and all (ti, . . . , t m ) in the support of /. Thus we can take the limit 
as N — > oo in J4.3| ) and obtain the right hand side of ( 1.14 ). This completes the 
proof. □ 

Before proving theorem 1.3 we need some preliminary results on the level spacing 
distribution. Let pn{x) be a symmetric probability density on WL N with correlation 
functions defined by (1.1). Assume that R^' /N — > p(t) (weakly) as N — > oo, so 
that p{t) is the asymptotic density. Let u be a given point such that p(u) > 0, 
and let t/y be a sequence such that t/y — > oo but t^/N — > as N — > oo. Set, for 
|r| < 1/2, 



RmH<Tl> ■ ■ -,^ m ;r) = 



1 



2tjvr + (Ti 



, M + 



2tj\rr + <7„ 



) 



(Np(u)) m " l v Np(u) ' iVp(u) 

and let lZ m (o~i, . . . , cr m ) be the limiting correlation functions, which we assume are 
continuous, symmetric and translation invariant. Assume that, for each s > 0, 



(4.4) 
Set 



D N (s) = ^— sup \TZ m (ax, ■■■ ,<J m ) \ < oo. 

L — ' 777, U.|<, 
m=N+l \ a i\^ s 



00 / -i \m p 

H ^ = 5Z M- / ^">i> • • • > 



'[0,s 



(the probability of no particle in [0, s]), which is well defined by (4.4). Also, set 



(4.5) 



AN) _ 



sup l^f^CTi, . . . ,<J m ;r) - 1l m (<Ji, . . . ,a m )\. 
kil<s,H<i/2 



Proposition 4.1. Let Sn(s,x) be defined by (l.l). Then 

< D N (s) + £ 



(4.6) 



SN(s,x)pN(x)d X 



H"(u)du 



N 



m — 1 



(m-1)! m ' 



Proof. We first show that 



(4.7) 



/" s W s m_1 /" 

/ H"{u)du = V] — - / 7?. m (0, r 2 , . . . , r m )dr 2 . . . dr„ 

Jo ^ 2 1)! yp),.]—! 
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see . Since lZ m is translation invariant and symmetric by assumption, we have 



(4.8) H>(u) = Vun-J2 

m 
oo 

= Km > 



i ^ 



(-l) m 1 



[-e,0] X [0,t»]* 



c— >0 z — ' TO,! £ 
m— 



V ( ~ 1)m / n (o 



(-^l ) • • • ; •Em)d % 



X2 ; ■ ■ ■ , %m)d X, 



where we have also used (4.4) and the continuity of 7?. m . Continuing in the same 
way we see that H(u) is act ually a C°° function, in particular H"(u) is well defined 
and continuous. From ( [4.8| ) we get 



(m- 1)! 

m=2 v y 



7?. m (0,a;2, . . . ,x m )d m 1 x. 



Hence H'(0) = —lZ m (0) and we see that the right hand side of (4.7) equals H'{u) 
-ff'(O), which is what we wanted to prove. 

It is proved in Q], using a result from |L7|, that 



Sn(s, x)pN{x)d N x 



/ dr 



-1/2 J[0,mm(s,(l-2r)tiv)] T 



Hence, the estimate (4.6) follows from (4.4), (Iha) and (4.7). 



□ 



We turn now to the proof of theorem 1.3. 



Proof. Just as in the proof of theorem 1.2 above we see that since P € W 6+e and 

N^jvlloo < JV/2*JV. 

(4.9) 



S N (s,x(M))dQ^(M) 



Ht 



S N (s,x)p N (x;y(H))d N x) dP^(H) 



< C—N 2 - {&+t)(1/2 -^ < — , 

if we take £ sufficiently small, and also that (B.4) holds. From proposition |1 . 1| , ( |2.4| ) 
and proposition |2.3| we know the correlation functions of pn{%', y), and if we take 
un = u + (2t^r + <Ji)(N )o(w)) _1 in lemma |3.2| we see that 

(4.10) 
1 



-IC(u 



2t N r + Ui 2t N r + a J 



sin7r(cri — Oj) 



Np(u) y Np{u) ' Np{u) 



<C(^+N-t) 
= u N 



for a. a. H [p( N )\. Thus, the limiting correlation functions are 

' sin7r(cr J ; — Oj) " 



Tl m (ai, . . .,a m ) - dct 



7r(CTi - <j,j) 



«.i=i 
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Since the matrix in the determinant is positive definite it follows from the Hadamard 
inequality that 



D N (s)< £ 

m=N+l 



Also, since 



K^\a 1 ,...,a m ;y) = det(-^ ) lC(u , 



2tNf + C; 2t]\rr + <7j 



-,u- 



Np{u) 



v) 



it follows from ( 4.1C ), the multilincarity of the determinant and Hadamard's in- 
equality that 

\n^(a;y)-TZ m (a)\ < m(l + ^ r^m" 1 / 2 , 



and hence em < m{l + lonY" 1 u>N'm 
formula and the fact that lun — > 0, 



»n/ 2 Now, by proposition 4.1, Stirling's 



(4.11) 



Sn(s, x)pN{x;y{H))d x — I H" (u)du 



< 



m=N+l 



V — + UJ N V 

z — ' m' z — ' 



iV 



(m- 1)! 



(l+Wjv) m_1 TO (m+2)/2 =0(1) 



as iV — > oo, for a. a. H [PW]. If we combine (gj) and we see that the 

theorem is proved. □ 
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